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ABSTRACT 

We consider a slab of nuclear matter and investigate the collective excitations, which 
develop in the response function of the system. We introduce a finite-range realistic 
interaction among the nucleons, which reproduces the full G-matrix by a linear com- 
bination of gaussian potentials in the various spin-isospin channels. We then analyze 
the collective modes of the slab in the S = T = 1 channel: for moderate momenta hard 
and soft zero-sound modes are found, which exhaust most of the excitation strength. 
At variance with the results obtained with a zero range force, new "massive" excita- 
tions are found for the vector-isovector channel . 
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1. Introduction. 

The slab of interacting nuclear matter has revealed itself as an interesting tool[l,2] 
for investigating specific properties of nuclei: being confined in one dimension only, it 
allows to consider, e.g., the influence of a surface (already present in the semi-infinite 
slab of Esbensen and Bertsch[3]) together with the effects of discrete levels, which are 
appropriate for a spatially confined system. But at the same time it keeps typical 
features of an infinite system, like for example density oscillations with phonon-like 
dispersion relation, together with translational invariance along the two unconfined 
directions. 

This model also offers a schematic (and rather extreme) situation to be compared 
with the one of highly deformed nuclei: as already pointed out in ref.[2], the observed 
splitting of the dipole resonance can be related with the analogous phenomenon which 
shows up in the collective modes (zero sound) of the slab. 

In this paper we consider, as in ref. [2] , the response function of a slab of interacting 
nucleons within the framework of the Random Phase Approximation (RPA): however, 
while in the previous work a schematic zero range force (of the Landau-Migdal type) 
was employed, here we utilize a more realistic nucleon-nucleon interaction, which is 
derived by suitably fitting previous G-matrix calculations [4,5]. The effective inter- 
action derived in the present work has the appropriate finite range in all spin-isospin 
channels and allows a more detailed investigation of the interplay between range of 
the interaction and size of the system. This point will be explored in connection with 
the collective modes of the slab, as derived within the RPA framework: indeed we 
will show how the corresponding response function is modified in the presence of finite 
range forces. 

In Section 2 we briefly describe the model, summarizing the appropriate formalism 
to deal with the response function of the slab; in Section 3 we give the derivation of 
the effective particle-hole interaction, which is considered in the various spin-isospin 
channels. Finally in Section 4 the resulting RPA response function and, specifically, 
the collective modes of the slab are presented and shortly discussed. 

2. Description of the model. 

We consider, as in ref. [2], a system of noninteracting fermions, confined in the 
domain < z < L by an infinite potential well; the single particle wave functions are 
then: 
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(2.1) 



where S is the area of the surfaces which set the boundaries of the system (the slab) , 
Xs and £t are two-component spinors in the spin and isospin space and the symbols 
|| and _L correspond to the parallel and perpendicular directions to the slab surfaces. 
In (2.1), since at the end we shall let S — > oo, fey is assumed to be a continuum wave 
number whereas k± is quantized according to 
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with n=l,2,3,..., the total particle wave number being k = \Jk 2 + k\. Eq.(2.2) also 
relates M with the slab thickness L and to the Fermi momentum fop: since M has to 
be an integer number, by fixing the slab thickness fop turns out to be quantized (or 
viceversa, indeed for convenience we will assume a fixed value for fcp = 1.36 fm -1 ). 

From the single particle wavefunctions and (kinetic energy) eigenvalues one can 
obtain the polarization propagator (or particle-hole Green function) for the non- 
interacting Fermi gas of nucleons in the slab; it reads 



n (xi, x 2 ) = -jG°(xi, x 2 )G°(x 2 , Xl ) 



(2.3) 



where G° is the free single particle Green's function and the factor of 4 arises from the 
traces over spin and isospin (we assume an isospin symmetric, N=Z, nuclear matter). 
The Fourier transform for the slab system is then defined as follows: 



f(k\\,kj_,k' ± ) = S / d r\\exp 
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f(r\\,z 1 ,z 2 ) = 
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x ^2 f(k\\,k±,k'_ L )cos(k_ L z 1 )cos(k / _ L z 2 )r]k ± Vk> ± (2.5) 
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where fey should be thought of as fey — fey since the system is translational invariant in 
the direction parallel to the surface and rjk ± = 1(1/2) for fej_ 7^ (k± = 0). By applying 
the above transformations to the polarization propagator (for details and explicit 
formulas see ref.[2]) we calculate the real and imaginary part of n^ ab (fc||, fej_, k' ±1 00). 

We remind that the imaginary part of the polarization propagator is directly 
related to the response function of the system, thus providing the excitation strength 
of the slab when this is coupled to an external probe; if the latter, for example, induces 
charge density fluctuations into the system, the response per particle is given by the 
relation: 

2 

R(k\\,k±,u) = -— ImU slab (k h k ± ,k ± ,u) . (2.6) 

As already stated in the introduction, in this work we are mainly concerned with 
the effects of a realistic interaction on the slab response function, and with the interplay 
between the finite range of the force and the finite size of the system. These items are 
best investigated within the framework of the Random Phase Approximation ( RPA 
) scheme. We have already pointed out in ref.[2] the main differences between the 
infinite nuclear matter and the slab RPA equations for the polarization propagator. 
Due to the finite size of the slab in the transverse dimension, translational invariance 
is lost and the simple algebraic RPA equation[6] for the polarization propagator in 
nuclear matter turns into a set of coupled equations for the non-diagonal elements of 

n^at^ * n momen t um space. 

In order to derive the latter it is convenient to start with the RPA equation for 
the polarization propagator in r-space; it reads (the subscript "slab" is omitted to 
simplify the notation): 



n RPA (fi - f 2 , z 1 , z 2 ) = n°(fi -f 2 ,z u z 2 ) + J d 2 r 3 d 2 r 4 J dz 3 dz 4 



x n°(fi - r 3 , z u z 3 )V(f 3 - f 4 , z 3 , z 4 )U RPA (f 4 - f 2 , z 4 , z 2 ) (2.7) 



where the r*j are two dimensional vectors in the (x,y)-plane. Since the interaction is 
time independent, we shall not explicitly indicate the time (or energy) dependence in 
our formulas. 

Applying now the cosine Fourier transformation, defined in (2.4) and (2.5), to 
both sides of (2.7) we get, after some manipulations, the compact expression 
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n (fe||,fe ± ,A; ± ) . (2.8) 



This equation is easy to solve numerically since the transverse momenta are dis- 
cretized according to k± = nkp/M (n = 0, 1, 2 . . .) and the longitudinal momentum 
fey can take arbitrary values and plays the role of a fixed parameter. For a definite 
choice of fey and cj (which is implicit in the free and RPA polarization propagators), 
eq. 2.8 becomes a two-dimensional system of linear algebraic equations. We cut the 
range of the allowed &x-values when it exceeds by more than three times the upper 
limit of the corresponding response region for the free system. Moving to k± values 
higher that does not affect the results by more than one percent and doesn't 

change the qualitative features of the solution. 

Concerning the particle-hole interaction V(k\\,k'j_,k , j_) between nucleons in the 
slab system, it is derived from a finite range realistic NN potential, as it will be ex- 
plained in the next Section. We anticipate here that in nuclear matter it is represented 
by a linear combination of gaussian functions both in momentum space, 



V{k)=ge- k2 ' m \ (2.9) 

where g and m are fitting parameters, and in r-space, where it reads, correspondingly, 

V(r) = Ce-^ r \ (2.10) 

with n = m 2 /4, C = g(m s /8n 2 )^n. 

For the slab system we have to perform the Fourier cosine transformation (2.4) 
on gaussian functions like (2.10) and the resulting particle-hole force (apart from the 
spin-isospin matrix elements) turns out to be of the following form: 

L L 

V{k\ b k'^k'[_)=SC 7 ^e k V^ J J d Zl dz 2 e-^ Zl - Z2)2 cos(k' ±Zl ) cos(k'iz 2 ) , (2.11) 







where the integration over z±, z 2 has to be carried out numerically. As a result of the 
finite range of the interaction V(k\\,k'j_,k'j_) is neither diagonal with respect to k' ± and 
k'[_ nor constant in momentum space (at variance with ref.[2]). 
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3. The Effective Interaction. 



To proceed with the calculation of the RPA response function, we need the finite 
range potential which acts between nucleons. We start from any realistic NN potential 
derived by fitting the NN scattering data and, in order to properly deal with the strong 
short range correlations, we insert this "bare" interaction in the Bethe-Goldstone 
equation. In the operator form the latter reads: 

G = V + V^ 1 -G (3.1) 

where V is the bare NN potential, G is the so called Brueckner G-matrix (which will 
coincide with our effective interaction), E is the starting energy, H is the Hamiltonian 
operator for the intermediate two-particle states; finally Q is the Pauli operator, which 
takes care of the medium effects of the Fermi system, by forbidding the particles to 
scatter into occupied intermediate states. 

The Bethe-Goldstone equation is then solved in the infinite nuclear matter, fol- 
lowing the method suggested by Haftel and Tabakin[7]. Eq.(3.1) implicitly depends 
upon the density of the system via the Pauli operator; we fix the starting energy to be 
E=74 MeV. Further details can be found in references [4,5]. As a solution we have the 
so called G-matrix, which can then be expressed in terms of the direct and exchange 
matrix elements of an effective potential: the latter is subsequently represented by 
means of a suitable parameterization. 

The method we employ has been developed in ref.[5]: working in the particle- 
hole representation, one interprets the resulting G-matrix as direct plus exchange 
matrix elements of some effective potential, whose dependence upon spin and isospin 
operators is explicitly taken into account. Then the momentum dependence of this 
potential is suitably parameterized in terms of simple Yukawa's functions, with the 
aim of obtaining an effective potential which is vaguely related to the exchange of 
some "effective" mesons. 

For the purpose of applications to the slab system the Yukawa's form of the in- 
teraction hinders the (partially) analytical evaluation of the particle-hole force (2.11). 
Thus in the present work we have considered a local parameterization of the G-matrix 
elements in the form of simple gaussians. These might turn quite advantageous also 
for finite nuclei calculations, e.g. when the single particle wave functions are expressed 
in terms of Harmonic Oscillator eigenstates. 

To illustrate the fitting procedure let us consider the particle-hole matrix elements 
of G in a definite spin-isospin channel, < ST\G(k,pi,p2)\ST >=< ST\U p h\ST >; as 
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shown in Fig.l, the direct matrix element depends on the total p-h momentum k while 
the exchange one depends upon pi — p 2 (the center of mass dependence being taken 
into account by some averaging procedure). 

We now assume that the momentum dependence of the G-matrix, solution of 
eq.(3.1) with some realistic input for the NN bare potential, can be accurately repro- 
duced by a two-body potential of the form: 

U ph (q) = goo e- q * lm *^o % + 9oi e"? • f 2 ) 

+ 010 e-« 2 / m ?o(^ ■ a 2 )l T + 9ll e^'™^^ ■ a 2 ){r 1 • f 2 ), (3.2) 

where for sake of simplicity we do not include tensor components, which could be 
extracted as well by considering definite spin projections in the G-matrix elements. 
Considering, for example, the 5 = 0, T = channel, the p-h matrix element of (3.2) 
reads (with the same notation of Fig.l): 

< m\U ph \m >= g 00 (4e" fc2/m o° - e~^~^ x 

- 3#oi e'^-v^K^ - Zg loe - { ^-v 2)2/m ™ - 9 gil e -^-V2) 2 lm\ K ^ 

With the approximation p\ — p 2 — (which turns out to be accurate, at least for not 
too large particle momenta) the exchange contributions simply reduce to a constant 
and, in order to fix the parameters in the effective potential U p h, eq.(3.2), one has to 
solve the following system: 

< 00|C(fc)|00 > = g 00 (4e" fc2 / m o° - l) - Sg 01 - 3g w - 9g n 

< 01|G(fc)|01 > = -g 00 + 9oi (4e- k2 / m 'n + _ 3gw + 3gu 

< 10|G(fc)|10 > = -goo ~ 3^oi + 9io (Ae~ k2 ' m ^ + l) + S 9ll 

< ll|G(fe)|ll > = -goo - 9oi + 9io + 9ii (le^^ - l) 

This allows to extract strengths (gsr) and masses (msr) for all components in U p h- 
In order to obtain a good fit to the original solution two gaussians are needed for each 
spin-isospin channel. Moreover we found it more convenient (and accurate) to find 
separate parameterization for the full G-matrix elements and for the direct part of 
them. For example, in the S = T = channel one sets: 

< 00|G|00 > tot = 4 (g&] F exp{-k 2 /(m^ F ) 2 } + gf \ F exp{-/ C 2 /(m| ) 2 ) )F ) 2 }) , (3.5a) 
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(3.4) 



< 00|G|00 > dir = 4 (gg] D eM-k 2 /(m£ ] D ) 2 }g$ D exp^/rf^) 2 }) . (3.56) 

The parameters of the present fit are summarized in Table I for the direct part 
of the G-matrix and in Table II for the full interaction. We observe that, at variance 
with ref.[5], the present parameterization involves the full G-matrix and not only the 
corrections induced by the ladder diagrams on the original bare potential (which in 
some cases are artificially large) ; thus one does not need the explicit knowledge of the 
NN interaction employed in the Bethe-Goldstone equation. The effective particle-hole 
potential resulting from our fit is illustrated in Fig. 2, where the full p-h interaction is 
displayed for each spin-isospin channel and compared with the exact G-matrix. The 
quality of the fit is fairly good, but for the scalar isoscalar channel, which however 
will not be utilized in the following: indeed in this case the interaction turns out to 
be strongly attractive over the whole range of momenta and cannot produce the hard 
collective modes we are interested in. 

4. Results and discussion. 

We have solved the RPA equations (2.8) for the slab polarization propagator by 
utilizing the effective potential of the previous Section. As one can see from Fig. 2 the 
full interaction develops a large repulsive p-h force in the spin-isospin channel (S = 
T = 1), together with a rather strong momentum dependence; the latter (although 
we are neglecting the tensor components of the interaction) should be ascribed to 
the presence, in this channel, of the long range one pion exchange potential. As it is 
well known [8] the attractive pion exchange modulates the strong repulsive short range 
correlations in the spin-isospin channel, thus producing some softening of the RPA 
response function. 

From the resulting momentum dependence of the p-h force we expect indeed to 
observe a different behaviour, in the collective response of the system, with respect to 
ref.[2], where a zero-range (constant in momentum space) interaction was used within 
the same theoretical framework. One should also remind that the interaction utilized 
here is not diagonal in the perpendicular momentum, thus involving a larger influence 
of the non-diagonal terms of IT in the RPA equations (2.8). 

Let us then focus on the vector-isovector channel: it is well known that, in nuclear 
matter, at small energy and momentum transfers the repulsive interaction is strong 
enough to produce a collective excitation (zero sound) outside the continuum p-h 
response, which in turn is severely depressed. 
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In the slab system and with the present finite range interaction we have found, 
in analogy with ref.[2], two distinct collective modes, according whether the transfer 
momentum is parallel to the infinite dimensions of the system or perpendicular to it: 
in Fig.3a,b we display these zero sound dispersion relations, which have approximately 
a linear form: uj = vk, v being the corresponding velocity of propagation of the bosonic 
excitation in the nuclear medium. Results are displayed for infinite nuclear matter 
and for slabs with two different thicknesses (which are fixed by the value of the integer 
M). Fig. 3b shows a magnified view of Fig. 3a at somewhat larger momentum transfers, 
where one can better distinguish the different modes occurring for the M = 2 and 
M = 8 slabs. 

More precisely, the sound velocity in nuclear matter turns out to be ij™ 771 - = 
1.62vf (it would be ij"- m - = 1.71vf with a 5-force of approximately the same strength 
as the one of the gaussian interaction at q = 0); vf = 0.286c is the Fermi velocity. 
For the M = 2 slab the two sound velocities obtained with our gaussian interaction 
are (in parenthesis the corresponding values for the 5-force) uy = 1.36^ (1.5vp) and 
v± = 1.54vp (1.89vp), for the parallel and perpendicular modes, respectively. In a 
thicker slab (M = 8) the resulting sound velocities are closer to each other and to the 
infinite nuclear matter value: v\\ = 1.56vf (1.Q7vf) and v± = 1.58vf (1.79vf)- 

With respect to the analogous outcome in ref. [2] one should notice the following 
qualitative differences : 

i) the dispersion relation in nuclear matter is modified by the present interaction, 
since as the momentum increases the p-h force weakens (see Fig. 2) thus producing 
some softening with respect to a constant interaction; 

ii) for a fixed thickness of the slab, one still finds two distinct collective modes, the 
faster of which being associated with a transverse momentum transfer ( "perpen- 
dicular" mode), however both the perpendicular and the parallel modes lie below 
the nuclear matter curve, while in the previous work they were sitting on oppo- 
site sides. This outcome should signal the role played by the off-diagonal (in the 
perpendicular momentum) terms of the interaction, which were absent for the 
5-force and produce a sizeable softening both of the longitudinal and, even more, 
of the transverse modes. 

iii) we remark that the above discussed softening is due to the attractive part of the 
interaction and this, in turn, has to be associated with pion exchange, namely 
with the longest range component of the NN force. It is thus worth noticing that 
the relative lowering of the transverse zero sound velocity is much larger than the 
one of the parallel mode: this fact signals an important correlation between the 
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range of the interaction and the transverse dimension of the slab, where the finite 
size of the system plays the major role. 

Beside the previously illustrated collective modes, the present model displays 
additional excitations, which in the RPA scheme are signalled by poles in ReU RPA , 
not occurring in IT , corresponding to very narrow peaks in ImIl RPA (and thus in the 
response function); this is shown in Fig. 4, for the M = 8 slab. The energy-momentum 
dispersion relation of these modes, which occur at higher energies with respect to the 
zero sound, differs from the latter since it does not vanish in the zero momentum 
limit and then stays almost constant, with a weak quadratic dependence. They are 
illustrated in Fig. 5, for the M = 2 and M = 8 slabs: in the thicker system the energy 
interval between any two of these modes is smaller (they appear to be almost equally 
spaced); however the strength in the corresponding peaks (see Fig. 4) rapidly goes 
down. Indeed they disappear in the infinite system. 

The origin of these collective excitations lies in the finite range of the interaction 
utilized here and critically depends upon the interplay between range of the interaction 
and extension of the system: they did not show up when a zero-range force was 
employed, as in ref. [2]. Moreover the peaks seem to be more pronounced when the 
range of the interaction (2-3 fm) is comparable with the transverse dimension of the 
slab: in the M = 8 slab (L ~ 18 fm) only the lowest energy mode displays a significant 
strength. 

This situation resembles the low frequency modes occurring in crystals and in 
spin systems (ferromagnets, antiferromagnets, etc.): in both cases one observes (one 
or more) acoustical and optical branches, the former having vanishing frequency when 
k — > 0, the latter requiring finite excitation energy in the long wavelength limit. In a 
ferromagnetic system the above mentioned excitations are spin waves usually called 
magnons[9]: experimental evidence for these modes has been revealed by inelastic 
neutron scattering[10]. From the microscopic point of view, they could be explained 
in terms of collective excitations of an infinite range Ising modelfll], which allows 
to associate the finite frequency of the optical branch to the order parameter of the 
ferromagnetic phase. 

In the present case we notice that the above mentioned "massive" modes have 
been found in the spin-isospin channel (S = T = 1) and may be explored by letting 
the system to interact with an external electromagnetic field. They arise only when 
the interaction has a non-zero range and disappear when translational invariance is 
restored, by letting L — > oo; indeed in infinite nuclear matter the zero sound alone is 
allowed as a collective mode, irrespective on the (finite or zero) range of the model 
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interaction utilized. It is worth reminding that the only physical example of a similar 
mode occurring in an infinite system, the electron gas, are the so-called plasmons, 
whose origin is driven by the infinite range of the Coulomb potential. Here instead 
we have finite range forces and the "plasmon-like" excitations survive only as far as 
the range of the interaction is comparable with the slab thickness L. Actually in 
the infinite system the three different types of collective excitations which we have 
found (parallel and perpendicular zero sound, and plasmon-like) collapse into a single 
collective mode, with linear dispersion relation. 

It might be interesting to notice that in a Wigner lattice the occurrence of a 
similar situation has been interpreted in connection with the phenomenon of spon- 
taneous symmetry breaking. This model displays three different types of elementary 
excitations: two transverse modes which behave as phonons, with uj oc q, and a longi- 
tudinal collective mode (the "plasmon") which has a finite frequency as q — > 0[12-14]. 
According to Anderson the plasmon can be viewed as the massive "Higgs" boson 
associated with the spontaneous breaking of gauge invariance which is due, in the 
presence of the Coulomb potential, to the density fluctuations in the electron gas. 
Notably, the existence of this massive mode is crucially related to the long range of 
the Coulomb interaction. Short range correlations, like in crystals, would give rise to 
massless Goldstone bosons alone. 

As a final remark we would like to point out that the present calculation, although 
rather crude in the description of the nuclear confinement, yet utilizes a realistic p-h 
force, thus allowing to bear some confidence in the qualitative (if not quantitative) 
results we obtain for the collective excitations of the system. In particular we believe 
that the analysis of the interplay between the range of the interaction and the finite 
size of the system might offer some hints also for experimental research in highly 
deformed nuclei. 
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Table I 



Direct Matrix Element 


S T 


9 ST 

[MeVfm 3 ] 


\V ) 

77% ST 

[MeV] 


9 ST 

[MeVfm 3 ] 


\2 ) 
771/ s'jp 

[MeV] 





-1014.40 




4379.0 


423.25 


313.6 


1 


187.45 




449.7 


425.23 


2773.5 


1 


-426.00 




4478.2 


163.00 


654.8 


1 1 


-266.24 




9553.0 


272.70 


370.8 



Table I - Parameterization of the Direct p-h matrix elements of the G-matrix of 
ref.[4,5]; S,T denote the total spin and isospin quantum numbers of the relevant 
channels. 

Table II 



Direct + Exchange Matrix Element 


S T 


Qst 
[MeVfm 3 ] 


77% Q rp 

[MeV] 


(2) 
9 ST 

[MeVfm 3 ] 


(2) 

77% Qijp 

[MeV] 





-2865.64 


862.0 


185.90 


862.0 


1 


576.20 


10000.0 


576.2 


377.0 


1 


321.14 


9691.0 


606.16 


669.8 


1 1 


1307.6 


461.0 


-388.74 


1363.0 



Table II - Parameterization of the Direct+Exchange p-h matrix elements of the G- 
matrix of ref.[4,5]; notations are the same as in Table I. 



Figure Captions 



Fig. 1 - Graphical representation of the direct (a) and exchange (b) particle-hole 
matrix element of the interaction. 

Fig. 2 - The present gaussian fit (continuous lines) to the p-h matrix elements of 
the full (direct plus exchange) G-matrix (dashed lines) is presented as a function 
of k (in fm -1 ) for the various spin-isospin channels. The units on the vertical 
scales are MeV fm 3 . 

Fig. 3 - The zero sound dispersion relations for the infinite nuclear matter (con- 
tinuous line) and for two slabs with M = 2 and M = 8: the parallel modes 
correspond to the short-dashed (M = 2) and dotted (M = 8) lines, while the 
transverse modes (which are not visible in upper part (a) of the figure) are rep- 
resented by the dashed and dot-dashed curves, respectively. The latter must be 
considered as interpolations, to guide the eye, of the discrete-g points where the 
transverse zero sound can be found. 

Fig. 4 - Reil RPA (dotted line) and Imil RPA (continuous line) for the M = 8 slab, 
at <7|| = 0.4 fm -1 , q±_ = 0, versus energy. 

Fig. 5 - Energy versus momentum behaviour of the collective modes in the parallel 
direction for the M = 2 slab (dashed lines) and the M = 8 slab (dotted lines); 
the infinite nuclear matter zero sound is also displayed (continuous line). 
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